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Abstract 

' Charmed mesons may be produced when a primary cosmic ray or the leading hadron 

' in an air shower collide with an atmospheric nucleon. At energies > 10 8 GeV their 

decay length becomes larger than 10 km, which implies that they tend to interact in 
the air instead of decaying. We study the collisions of long-lived charmed hadrons 
in the atmosphere. We show that (A c , £))-proton diffractive processes and partonic 
collisions of any q 2 where the charm quark is an spectator have lower inelasticity 
than (p, 7r)-proton collisions. In particular, we find that a D meson deposits in each 
interaction just around 55% of the energy deposited by a pion. On the other hand, 
collisions involving the valence c quark (its annihilation with a sea c quark in the target 
or c-quark exchange in the t channel) may deposit most of the D meson energy, but 
their frequency is low (below 0.1% of inelastic interactions). As a consequence, very 
energetic charmed hadrons may keep a significant fraction of their initial energy after 
several hadronic interactions, reaching much deeper in the atmosphere than pions or 
protons of similar energy. 



1 Introduction 



Cosmic rays reach the Earth with energies of up to 10 11 GeV. When they enter the atmo- 
sphere they experience ultrahigh energy collisions, probing a scale not accessible at colliders. 
These processes are interesting because they could produce exotic particles or involve inter- 
actions not seen at lower energies. There is, however, a second generic reason that make these 
cosmic-ray collisions interesting: they probe a regime where the properties of the standard 
particles may be substantially different. 

In particular, here we will focus on the charm quark. Charm-quark pairs may be abundant 
inside extensive air showers [HE]- Once produced, each quark will result into a D +, ° meson 
or a A+ baryon, that will then decay weakly into light hadrons and (with a ~ 0.1 branching 
ratio) leptons. These processes have been extensively studied in the literature, as they may be 
the dominant source of atmospheric muons and neutrinos of energy E > 10 6 GeV [31 HJ El E] • 
A new regime, however, would be achieved if the charmed hadron has an energy E > E c 
such that its decay length (A^ = crE/m) becomes larger than its interaction length (Aj) in 
the atmosphere. For example, a 10 s GeV D + meson is effectively a long-lived (A^ + ~ 18 
km) heavy hadron that interacts and propagates in the air instead of decaying. Such a 
behaviour, unknown of in colliders (where a D meson decays before it has experienced any 
type of interaction) could look somewhat similar to the one of the i?-hadrons predicted in 
some supersymmetric models [3 [8]. The propagation in matter of a charmed hadron is of no 
interest at the Tevatron or the LHC, but it is most likely important in extensive air-shower 
experiments. 

In this article we analyze the collisions of very energetic D mesons and A+ baryons with 
protons at rest. Our objective is to estimate the amount of energy that they deposit in 
each hadronic interaction in the atmosphere. First we study diffractive processes, using the 
Montecarlo code PYTHIA [9] to simulate the collisions. Then we consider partonic collisions 
of any q 2 , with our analysis also based on PYTHIA. Finally we estimate the effects of going 
from a proton to a nucleus target. 

2 Diffractive processes 

A charmed hadron H is basically different from a pion or a proton in the fact that it 
contains a heavy core of mass m c ~ 1.27 GeV. If one thinks of a proton as three clouds of 
mass m = 0.3 GeV associated to the three constituent quarks, then the A c baryon consists 
also of three similar clouds, but one of them has an additional electroweak core and a total 
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mass m c + m — 1.6 GeV. In a long-distance (low-g 2 ) collision with an atmospheric nucleon, 
however, this heavy core will be invisible (the proton and the A c clouds will look identical), 
as only interactions of q 2 > m 2 c can resolve it. Therefore, we will assume that the momentum 
exchanged through pomerons or other non-perturbative dynamics with the target does not 
depend on the electroweak core in H. 

Let us be more specific. We separate the mass m# of the heavy hadron into 

m H = m c + A . (1) 

The light degrees of freedom in H carry just a fraction w = (ran — Tn c )/mH of the hadron 
energy E, so in a diffractive process H will be seen by the target nucleon like a light hadron 
of energy wE. Therefore, to estimate the momentum g M absorbed by H in the process, we 
will just simulate with PYTHIA the collision of a proton (for H = A c ) or a pion (for H = D) 
of energy wE with the nucleon and will assume that is the same when the incident particle 
is the charmed hadron. 

Once H absorbs the momentum g M it becomes a diffractive system H&f of mass M, the 
critical parameter in the collision. If M < niu + 1 GeV the collision is quasielastic and H^f 
will just decay into two bodies {e.g., H + if). For larger values of M the system is treated 
by PYTHIA pU] like a string with the quantum numbers of H . When if is a baryon the 
string may be stretched between a quark and a diquark or between a quark, a gluon and a 
diquark, whereas for a diffractive meson the string connects a quark and an antiquark or a 
quark, a gluon and an antiquark. 

To illustrate our procedure, let us consider the diffractive collision of a 10 9 GeV A+ 
baryon (m^ = 2.28 GeV) with a proton (m p = 0.94 GeV) at rest. We need to simulate 
the collision of a 4.43 x 10 8 GeV proton, and we will do it in the c.o.m. frame, where each 
proton carries 14.4 TeV. 

(i) In a first PYTHIA example the diffracted proton absorbs g M = (—1.3; —0.13, 0.10, —1.3) 
GeV, getting a mass of M — 1.58 GeV. The system then decays into a A + 7r°, and the 
A + finally produces a proton and another 7r°. Going back to the lab frame we find the 
final proton with E' = 2.75 x 10 8 GeV, i.e., the leading baryon carries a fraction 0.62 
of the initial energy. The elasticity z in this collision is then z = 0.62. If the incident 
particle is a 10 9 GeV A+, we first go to the same pp c.o.m. frame (where the energies 
of the charmed baryon and the target proton are 32.5 TeV and 14.4 TeV, respectively). 
The A+ absorbs there the same momentum q^ and gets a diffractive mass of M = 3.2 
GeV, it then goes into S+7r°, and the E c decays into a A+ and another n°. Going 
back to the lab frame we obtain that the final A+ baryon carries E' = 9.3 x 10 8 GeV, 
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implying an elasticity z = 0.93. This type of quasielastic processes accounts for 8% of 
all diffractive collisions. 

(ii) In a second PYTHIA simulation the momentum absorbed by the incident proton in 
the c.o.m. frame is = (0.12; 0.14, 0.14, —0.12) GeV, which produces a diffractive 
mass M = 83 GeV. The system evolves into a string stretching between a w-quark and 
a (ud)-diquark. Quark fragmentation and the (strong or electromagnetic) decay of 
baryonic resonances results then into a leading baryon plus 18 other hadrons of lower 
energies. In particular, in the lab frame there is a p of 1.8 x 10 8 GeV (i. e.,z — 0.41) plus 
16 mesons (pions and kaons) and a, pn pair sharing the rest of the energy. Changing 
the proton for a 10 9 GeV A+ we obtain a diffractive mass of M — 124 GeV. Now the 
system may define 3 different (equally probable) diquark-quark strings: (cd)-u; (cu)-d 
or {ud)-c. The first two cases tend to result into a leading charmed baryon after the 
collision. For example, in the first case we obtain a final A+ of 6.0 x 10 8 GeV (z = 0.60). 
The third case, however, is basically different, as the struck quark is the charm and 
it will most likely fragment into a D meson. With the PYTHIA simulation we obtain 
a final D + of 3.8 x 10 s GeV (z = 0.38). In around 30% of diffractive collisions the 
A c changes into a D meson. In contrast, in a diffractive proton-proton collision the 
leading hadron becomes a meson in just 15% of the cases. 

(Hi) In the final example PYTHIA provides an event with large diffractive mass (M = 6087 
GeV; q» = (641; 0.36, 0.43, -641) GeV) and a diquark-gluon-quark string. The final 
spectrum is similar to the one in case (ii), although with a larger multiplicity of final 
states and a lower energy in the leading hadron. Whereas in the pp collision we find a 
neutron carrying a 38% of the initial energy, in the A c p process we obtain a 4.1 x 10 8 
GeV (z = 0.41) A+ or a 7.3 x 10 s GeV (z = 0.73) D° depending on the flavor of the 
struck quark. 

In Fig. 1 we plot the distribution of the fraction z of energy taken by the leading baryon 
after a pp diffractive collision or by the charmed hadron after the analogous A c p process. 
The average values of the distributions are (z) = 0.60 and (z) = 0.77, respectively, which 
imply an inelasticity K — 1 — (z): 

7^ = 0.40; tff p = 0.23. (2) 

The regions of highest z (with = 0.01) correspond to 28% of diffractive processes 

where the incident hadron does not break (but the target proton does), whereas 8% of 
processes correspond to quasielastic events with M < m^ c + 1 GeV and have a slightly higher 
inelasticity, K { £ p = 0.09 ± 0.04, where AK = J(z 2 ) - (z) 2 . In the remaining collisions the 
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Figure 1: Fraction z of energy taken by the leading baryon or the charmed hadron in pp and 
A c p diffractive collisions, respectively We have separated the distributions in the 4 types of 
A c p processes described in the text. 

A c may become a D meson or a charmed baryon (29% and 35% of all diffractive processes, 
respectively). The D mesons carry a fraction (z) = 0.64 of the initial energy (i.e., K\L = 
0.36 ±0.16), whereas charmed baryons imply a slightly lower inelasticity, K^ p = 0.32 ±0.17. 
The leading D meson may be a D + , a D° or a D s in an approximate proportion of 1 : 3 : 0.5 
(the D + -D° isospin symmetry is broken by the decay of charm resonances), whereas the 
baryons are mostly A+ (94%) with some S° (3%) and (3%). In Fig. 1 we have separated 
these four types of diffractive A c p processes. 

For an incident 10 9 GeV D + meson we proceed in an analogous way, studying the diffrac- 
tive collision of a 3.2 x 10 8 GeV pion with a proton at rest. We simulate with PYTHIA 
the collision in the c.o.m. frame, we read the absorbed by the pion and assume that the 
momentum absorbed by a D + in that frame would be the same. We find that in 24% of 
the processes the incident meson remains unbroken (i.e., M = mo), whereas quasielastic 
events (M < mp ± 1 GeV) account for 15% of the total. The inelasticity of these two types 
of events is K^J = 0.02 and = 0.05 ± 0.05, respectively. The rest (61%) of diffractive 
collisions deposit an average fraction Kj® = 0.33 ± 0.17 of the initial energy. In Fig. 2 we 
give the z distribution for the three types of events. The peaks in D collisions appear at 
z ~ mr,/(mr, + m ViU! , ...), whereas the step in the pion distribution results from processes 
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Figure 2: Fraction z of energy taken by the leading baryon or the charmed hadron in 717? and 
Dp diffractive collisions, respectively We have separated the distributions in the 3 types of 
Dp processes described in the text. 

with two final pions of similar energy. The average inelasticity is 

Kfl = 0.47 ; K% = 0.21 . (3) 

Whereas a pion loses energy faster than a proton, the diffractive collisions of a D-meson are 
slightly less inelastic than the ones of a charmed baryon. 

Finally, the diffractive cross section for a 10 9 GeV charmed hadron can be read from the 
PYTHIA estimate for a 4.4 x 10 8 GeV proton or a 3.2 x 10 8 GeV pion, 

4f = 26.2mb; 0% = 15.6 mb . (4) 



3 Partonic collisions 

Non-diffractive processes dominate the inelastic cross section in PYTHIA simulations. In 
particular, at 4.43 x 10 8 GeV the cross section for parton-parton interactions of any q 2 in pp 
and up collisions i£] 

^ = 62.2 mb ; <^ = 43.6 mb . (5) 
1 PYTHIA extrapolates to q 2 < 2 GeV using a minimum-biased method. 
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We will use these simulations to obtain an approximate description of partonic Hp collisions. 

We model H as a system with the same parton content as the corresponding proton 
or pion but substituting a valence up quark u for the charm quark. Like in diffractive 
processes, we will associate a charmed hadron H of energy E to a light hadron of energy 
E(rriH — m c )/m,H- If «o carries a fraction x of the proton or pion momentum we will change 
it for a c with 

m c rriH — m c 

x c = 1 x. (6) 

m H m H 

In this way the excess of energy in H is carried entirely by the charm quark, whereas the 
light partons in both hadrons (H and p or n) carry exactly the same amount of energy. 

We will then distinguish two types of Hp non-diffractive collisions: those where the 
charm is an spectator (i.e., it is a light parton in H who hits a parton in the target proton), 
and processes where the charm itself interacts. For the first one we will just simulate with 
PYTHIA the parton process using a light hadron and then substitute the spectator uq for 
the charm quark. Charm interactions, on the other hand, have a much smaller cross section, 

& = 0.8mb, (7) 

than the processes with an spectator charm quark, 

a c A ^ ec = 61.4mb; a% ec = 42.8 mb , (8) 

but they imply collisions of higher inelasticity. In particular, there is the possibility that the 
c quark in the incident H hadron annihilates with a sea c in the target proton. The charmed 
hadron after such process has lost basically all its energy. These events, however, occur in 
just a fraction of all partonic processes, 

ag p = 0.02 mb . (9) 



Let us first describe A c p non-diffractive collisions. In Fig. 3 we plot the fraction z of 
energy carried by the charmed hadron after the collision. We obtain that the leading hadron 
may be a baryon (55%) or a D meson (45%). The baryon may be a A+ or, with much smaller 
frequencies, or DP C . The three meson species [D + , D°, D s ) appear with approximate 

frequencies of (1 : 3 : 0.5). In Fig. 3 we have separated the spectra in the two cases. The 
average fraction of energy taken by the baryons is (z) = 0.60 (i.e., K ( ^ p = 0.40 ± 0.17), 
whereas for the mesons it is just (z) = 0.47 (K\J = 0.53 ± 0.15). The events where the 
charm interacts with a light parton in the target proton are mostly included among the ones 
with a leading D meson, and have an average elasticity of (z) = 0.42. Finally, in 0.03% of 
the partonic collisions the incident charm anihilates (it is traded by a spectator sea charm 
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Figure 3: Fraction z of energy taken by the leading baryon or the charmed hadron in pp 
and A c p non-diffractive collisions, respectively. We have separated the distributions when 
the final state is a charmed baryon (5) or a charmed meson (6). 

in the target proton) and K A ^ p m 1. The average inelasticity of non-diffractive collisions is 
(we include for comparison the inelasticity in pp collisions) 

K?- dif = 0.66 ; Kl~ dif = 0.46 . (10) 



Dp non-diffractive collisions are summarized in Fig. 4. The final charm hadron is almost 
always (85% of the events) a D meson carrying a fraction (z) = 0.55 of the initial energy 
(K^ p = 0.45 ± 0.16). The c quark anihilates (Jfjg » 1) with a sea c in the target just 
in 0.04% of the Dp partonic collisions. The inelasticity in these non-diffractive processes is 
substantially lower than in pion collisions, 

K^~ dif = 0.77 ; K n D : p djS = 0.45 . (11) 



4 //—nucleus collisions 



From the analysis in previous sections it results that the inelastic cross section for the collision 
of a 10 9 GeV charmed hadron with a proton at rest is 



a 



AcP = 88.4 mb ; a Dp = 59.2 mb , (12) 
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Figure 4: Fraction z of energy taken by the leading baryon or the charmed hadron in 7rp and 
Dp non-diffractive collisions, respectively. 




Figure 5: Fraction z of energy taken by the leading baryon in pp collisions, the leading pion 
in 7rp collisions, and the charmed hadron in A c p and Dp collisions. The distributions include 
diffractive and partonic collisions. 
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where diffractive processes contribute a 30% in A c p collisions and a 26% in Dp interactions. 
The cross section for proton and pion collisions of the same energy is around a 7% larger, 

a pp = 94.1 mb ; a np = 64.8 mb . (13) 

The average inelasticity in these collisions is then 

K AcP = 0.39; K Dp = 0.38, (14) 

which is substantially lower than the one in proton and pion collisions, 

K pp = 0.59; K np = Q.70. (15) 

These results depend basically on the fragmentation model used by PYTHIA, so they should 
not change substantially if the analysis were based on the Montecarlo code SIBYLL pj] or 
on any other code using the Lund string scheme. 

In this section we would like to comment on the approximate effects of going from a 
proton to a nucleus target. The total cross section for the collisions of H with an atomic 
nucleus of mass number A can be estimated as 

a HA ^ A 2/3 a H P ^ (16) 

where the factor of A 2 / 3 takes into account the screening between the nucleons inside the 
nucleus. For an averaged atmospheric nucleus of A = 14.6 this implies 

a A c air = 528 mb ; a D air = 354 mb . (1 7) 

The associated interaction length \f nt = m air / a H atr in the atmosphere when if is a charmed 
meson or a baryon is therefore 

A^ = 46g/cm 2 ; Af nt = 69 g/cm 2 , (18) 

which is a 6% and a 9% longer than those of a pion and a proton of the same energy, 
respectively. 

To deduce the spectrum of the leading hadron after the collision one must take into 
account that there may be more than one nucleon involved in each if -nucleus interaction. 
In general, it will be possible to distinguish between peripheral and central collisions. The 
first type is similar to the Hp processes discussed before, and we will then assume that the 
spectra coincide. Central collisions, on the other hand, imply a softer spectrum of secondaries 
and a smaller value of the fraction z of energy taken by the leading charmed hadron. 
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As a first order estimate, we will assume that peripheral and central collisions occur with 
equal frequency, and that the average inelasticity in central processes is the typical in a 
non-diffractive collision increased by 10%: 

K Hair « \ K%2r + I KSt « Q K Hp + i 1.1 Kg*) . (19) 
At E = 10 9 GeV we obtain 

K AcClir ^0A5; K Dair ^0AA. (20) 

For proton and pion collisions the same prescription gives an inelasticity 

K pair ^0M; K vaAr w 0.78 . (21) 

This 12% increase in K when going from a proton to a nucleus target compares well with 
the results obtained by other authors [12J. The frequency of baryon to meson transitions in 
A c -air central collisions would be close to the 45% obtained in partonic processes. 

Finally, it is important to estimate the probability p^ 5 air of a cc interaction that deposits 
all of the charmed hadron energy. Short-distance interactions in nucleus collisions scale 
proportional to the mass number, implying 

Aa Hp 

Hair S±°cc ,yy> 

Pec ~ A 2/3 a Hp ■ 

We obtain 

p A_ c air w Q.0015 ; pg air w 0.0020 . (23) 



5 Summary and discussion 

Charmed hadrons that decay weakly become long lived in the atmosphere at energies E > 10 8 
GeV, just like pions do at E > 10 2 GeV. The behaviour of pions in a calorimeter is well 
understood and, as a consequence, their dynamics inside extensive air showers is also known. 
We know, for example, that ultrahigh energy pions are produced in these showers but never 
reach the ground, since the atmosphere is equivalent to 10 meters of water vertically and 30 
times thicker from horizontal zenith angles. In contrast, D mesons are not directly observable 
at calorimeters in colliders because they decay before they can reach them. Any estimate of 
their possible effects in air showers must then first model its hadronic interactions. 

Here we have studied the collisions of charmed hadrons with protons. Our intention has 
not been to perform an analysis based on first principles, but to obtain an estimate of the 
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qualitative features of such processes based on PYTHIA. The charm quark inside the hadron 
carries a large fraction (m c /rriD ~ 0.7) of energy, so one may expect that the spectrum of 
the final state will be substantially different from the one in pion-proton collisions. We have 
argued that in most of the processes either the charm can not be resolved (q 2 <C m 2 c ) or it is 
just an spectator. This allowed us to correlate them with pion or proton collisions, and use 
then PYTHIA to simulate quark fragmentation and the (strong or electromagnetic) decay 
of higher-mass resonances. 

We obtain that the average inelasticity in Dp collisions is K^ p = 0.38, a factor of 0.55 
smaller than in irp interactions. If the target is an air nucleus, we estimate values around 
Koair = 0.44 and K^air = 0.78. These central values would imply that starting with an 
initial 10 9 GeV pion, after 10 hadronic interactions in the atmosphere one is left with leading 
meson of just 300 GeV. Instead, the energy of a D meson would be reduced to 3 x 10 6 GeV. 

In contrast to pions and protons, we find that the collisions of charmed mesons and 
baryons have a very similar inelasticity. However, A c baryons have a large probability of 
becoming D mesons after a couple of interactions. This may actually be an important 
source of D mesons in extensive showers if the intrinsic charm [I] dominates the PDFs at 
high energy in proton collisions. 

The proper inclusion of charm in extensive air-shower simulators should describe both its 
production (already in DPMJET [2]) and also its propagation. D and A c hadronic collisions 
are obviously irrelevant in colliders, which explains why they are absent in codes like SIBYLL 
pTj or QGSJET [13]. We think, however, that the possibility to search for observable effects 
in astroparticle experiments should make worthy an effort in such direction. 
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